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The results of numerical experiments on the structure of chaotic attrac- 
tors in the Khalatnikov - Kroyter model of two freedoms are presented. This 
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model was developed for a qualitative description of the wave turbulence of 
the second sound in helium. The attractor dimension, size, and the maximal 
Lyapunov exponent in dependence on the single dimensionless parameter F 
Q ' of the model are found and discussed. The principal parameter F is similar 

■ to the Reynolds number in hydrodynamic turbulence. We were able to dis- 

cern four different attractors characterized by a specific critical value of the 
parameter (F = F cr ), such that the attractor exists for F > F cr only. A 
simple empirical relation for this dependence on the argument (F — F cr ) is 
presented which turns out to be universal for different attractors with respect 
to the dimension and dimensionless Lyapunov exponents. Yet, it differs as 
to the size of attractor. In the main region of our studies the dependence of 
, all dimensionless characteristics of the chaotic attractor on parameter F is 

very slow (logarithmic) which is qualitatively different as compared to that of 
a multi-freedom attractor, e.g., in hydrodynamic turbulence (a power law). 
<^ . However, at very large F ~ 10 7 the transition to a power-law dependence has 

been finally found, similar to the multi-freedom attractor. Some unsolved 
problems and open questions are also discussed. 

1 Introduction 

x. 

In the present paper we continue numerical experiments ]I| with the Khalatnikov - 
Kroyter model |2| developed for a qualitative description of the wave turbulence of 
second sound in helium. 

The model is specified by an effective non-Hermitian Hamiltonian: 

H(ai, a 2 ) = {oJ x - i7i)|ai| 2 + (u 2 - 112)^2? + ^aja* 2 + /a* + c.c.) (1.1) 

where we slightly changed the notations in Ref.|2|, |IJ. This Hamiltonian describes 
the two linear oscillators via complex phase-space variables aj (j = 1,2) and the 
frequencies Uj — ijj with phenomeno logical dissipation parameters jj. Two other 
parameters of the model represent a nonlinear coupling fi of the two oscillators, and 
the driving force /. 
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The motion equations 

ihi = (ui - iyi)ai + 2fia{a 2 + f ^ ^ 

ia 2 = (co>2 - ^72)^2 + na\ 

were numerically integrated together with the corresponding linearized equations 
(for details see Ref.jJTJ). Particularly, all four Lyapunov exponents (Ai > A 2 > A 3 > 
A4) were computed. One of them, whose eigenvector goes along the trajectory, is 
always zero while the sum of all 

4 

Ar = ^ n = _2 (7i + 72) = const (1.3) 

n=l 

is the constant rate of the phase space volume contraction. 

Surprisingly, at a relatively weak force / this most simple model does describe 
the birth of turbulence (dynamical chaos) in a physical system. This was the main 
subject of studies in Ref.[]2j, [[J. 

As the force grows the two-freedom model is losing any relation to the real 
physical turbulence which is a multi(infinitely)-freedom phenomenon. Nevertheless, 
it seems of some interest, for the general theory of dynamical systems, to compare the 
structure and properties of a chaotic attractor in such an opposite limit represented 
by a sample model under consideration. It was the main motivation for us to 
continue general studies of the model on unrestricted range of its parameters. 



2 Scaling 

In the spirit of the theory of turbulence we introduce, first, a unique dimensionless 
parameter, similar to the Reynolds number, which determines all the dimensionless 
characteristics of the motion. To this end, we choose /, //, 7 as the three basic 
parameters, and form the desired combination 

F = t± (2.1) 

T 

which has the meaning of a dimensionless force. For this basic characteristic of the 
model were unique we need to fix all other dimensionless parameters: 

— = 0, — = 0.25, and — = -12.5 (2.2) 

^2 72 72 

The values of all three ratios are chosen as in the previous studies |2|, Particularly, 
we keep u\ =0, as well as the motion initial conditions 

ai(0) = a 2 (0) = (2.3) 
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This allows us to directly compare our new results with the former ones. In what 
follows we set 7 = 72 in the basic relation (2.1). 

Invariance of F with respect to variation of the three basic parameters gives us 
an extra freedom for choosing the latter in such a way to minimize the computation 
errors, and thus to reach higher values of F. 

In the present studies we were primarily interested in the properties of chaotic 
attractor. One of its principal characteristic is the metric (M) that is the set of 
signs of the Lyapunov exponents. As in the previous studies we did observe only 
one metric of the three possible, namely 

M = (+, 0, -, -) (2.4) 

with a single zero exponent A2 = 0. The next, more interesting, characteristic 
of attractor is the fractal (noninteger) dimension. By now, there is a dozen of 
various definitions for such a dimension (see, e.g., Ref.||). From physical point of 
view those all are meaningful and acceptable in principle. However, in numerical 
experiments the two of them are much more preferable. These make use of the 
Lyapunov exponents only which greatly simplifies the computation. We have chosen 
one, due to Kaplan and Yorke, because it is more close to other definitions [|, ||. It 
is given by the relation 

d = m + ^ n=1 , (2.5a) 
|A m +i| 

where m is the largest integer for which the sum J2™=i — 0. This is the simplest 
and widely used method for calculating dimension. In the model under consideration 
the dimension of a chaotic attractor is always within the interval (2 < d < 4). The 
upper bound corresponds to the dimension of the whole phase space, while the lower 
is the condition for chaos (A x > 0, see Eq.(2.4)). 

The second characteristic of a chaotic attractor, closely related to (but different 
from) the former is the maximal Lyapunov exponent Ai. In dimensionless form it is 

Ax = — (2.56) 

7 

Accordingly, the dimensionless sum of all the exponents Ar = Ar/7 = —2.5 (see 
Eq.(1.3)). " 

Besides, we computed two other, geometrical, characteristics of attractor: the 
average size and shape. The former is represented by attractor's rms radius in phase 
space 

R = y/Ri + R\ , R) = <A)> - <A j > 2 ^< A) > (2.5c) 

where Aj = djJ fJ>/ f (j = 1,2) are dimensionless variables, and the brackets <> 
denote the time average along a trajectory. In all cases the average shift of attractors 
< Aj > was relatively small. The attractor shape is characterized by the ratio 

s = r\ {2M) 
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The results are presented and discussed in the next Section. 



3 Results and discussion 

In Fig.l the attractor dimension d, Eq.(2.5a), is plotted in dependence on the prin- 
cipal parameter F, Eq.(2.1). The huge range of F comprises almost five orders of 
magnitude, from the first chaos border at F — F cr = 173 up to F ps 10 7 limitted 
by the computation time (about one full day). Five separated groups of points are 
clearly seen followed by a number of scattered points. Each group corresponds to 
a chaos window in F separated from the neighbours by the limit-cycle windows. 
Apparently, the latter also exist in the region of scattered points which all belong 
to chaotic attractors. However, these points are too few, because of computation 
difficulties, for the reliable location of any windows beyond F ~ 10 4 . 




ln(l+F) 



Figure 1: Chaotic attractor dimension vs. principal parameter F in a semilog scale 
(see text). Gaps between groups of points correspond to limit-cycle attractors. 
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The first two groups, studied already in our previous work [|l|], certainly belong to 
different attractors with different chaos borders F cr »s 173 and 190, respectively. It 
is natural, again in the spirit of the theory of turbulence, to introduce the principal 
argument of all dependences in the form: 

5F = F — F cr (3.1) 

where F cr is the corresponding chaos border which is characteristic, and different, 
for a particular attractor. 




Figure 2: Dependence of the attractor dimension d (upper points), maximal Lya- 
punov exponent A x (crosses), and — A 3 (lower points) on the parameter SF (3.1). In 
this figure F cr = 173 and 190 for the two first groups of points in Fig.l, and F cr = 
for all other points (cf. Figs.l and 3). Two solid horizontal lines indicate the di- 
mension interval (2 < d < 4). The dashed line is the mean < — A 3 >^ 1.25. Two 
curves are empirical relations (3.3) and (3.2) for d and Ai, respectively; F = 8.5. 

The result of such a rescaling is shown in Fig. 2: both groups have overlapped 
following some joint dependence d(SF). 
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The first surprise was in that one of the Lyapunov exponents proved to be nearly 
constant in the whole range of the empirical data: A3 « —1.25 ~ Ar/2. Even 
though the origin of this peculiarity is not known it allows for a simple relation of 
both empirical dependences in Fig.2. Besides, we did use this result as an additional 
check of the computation accuracy. 

Indeed, using Eqs.(1.3), (2.4), and (2.5a) we obtain a simple relation 

J-A 3 (d-2), d<3 
Al " 1 -A 3 -A r f5i d > 3 (3 - 2) 

with constant parameters A3 « —1.25 and Ap = —2.5 ~ 2A3 whatever the depen- 
dence d(SF). 

The next, and the main, result was a very slow dependence on F of all the 
dimensionless characteristics of chaotic attractors (see Figs.l and 2). This is why 
we make use of a logarithmic argument Fl = In (1 + SF) rather than of SF itself. 
That slow dependence prevents us from reaching a close vicinity of the upper bound 
d = 4, and thus from clarifying an interesting question in Ref.JI]]: if the dimension 
does asymptotically approach this bound or, instead, saturates at some lower value 
of d (for farther discussion see next Section) . The maximal value actually reached in 
our numerical experiments was approximately d ~ 3.73 only due to a rapid increase 
of the computation time as d — > 4. 

The empirical dependence d{F£) (and hence Ai(Fl) as well) has been constructed 
using the data for the two first groups in Fig.l with known F cr values. To this end, 
let us assume, first, that d(F L ) does reach the upper bound as F L — ► 00. Further, 
Fig.2 suggests that this dependence is close to linear near F L = 0, hence we added 
the unity under logarithm in Fl. Finally, we use a simple relation: 

d(F L ) = 2+ A f L F L = ln(l + 5F) (3.3) 

with a single fitting parameter Fq. 

Surprisingly, the simple empirical relation (3.3) fits, with a particular value of 
the fitting parameter F = 8.5, not only the points of the first two groups but also, 
with a reasonable accuracy (especially for d), the rest of points with unknown F cr 
values set to zero in Fig.2. Notice that deviations of empirical points from a smooth 
dependence are not due to computation errors, which were carefully checked, but 
apparently represent some fine structure of chaotic attractors we did not study in 
the present work (cf. Ref. []]]). 

Still, the agreement between empirical data and analytical relation is not yet 
satisfactory especially for the three rightmost points with the biggest F. 

For improving the agreement we used two methods. First, we tried to fit the 
unknown F cr values for the three groups of points around F L = 7 (see Fig.l and 
the Table). Interestingly, the two groups seem to belong to the same attractor even 
though they are separated by a wide gap with limit cycles only. Apparently, this is 
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because of the fixed initial conditions (2.3) which may or may not find themselves 
within the attractor basin for a particular value of Fl (see also Ref.[0, |1|]). For this 
reason we cannot directly see the vicinity of chaos border for these three groups of 
points. Hence, the two new values of F cr are the fitting parameters unlike the first 
two. 

Table. Characteristics of few-freedom chaotic attractors 



attractor's 


F 

cr 


F cr jF 


R(F) 


S(F) 


number 


chaos border 




size 


shape 


1 


173 


1 


1.41 


0.74 


2 


190 


1 


1.43 


0.75 


2 


190 


0.32 


1.7 


0.9 


3 


1913 


0.83 


1.9 


1.0 


4 


3038 


0.93 


2.0 


1.1 


F > 4000 


0.8 F 


0.8 


2-6.5 


1 - 1.6 



Beyond the five groups we managed to compute just a few scattered points which 
prevented the location of other attractors. In this region we used another method 
based on a simple assumption that each point belongs to a separate attractor whose 
chaos border is directly related to the position of this point 



F cr (F) = F s -F (3.4) 

where F s is a new fitting parameter. This method of a statistical nature is, in a 
sense, opposite to the first one. The assumption (3.4) is very crude, of course, but it 
allows us to farther improve the agreement with empirical relations (3.2) and (3.3) 
including two points at Fl ~ 13.5 but not the last one at the biggest Fl « 14.5. 

Certainly, a couple of these 'abnormal' points is too few for any definite conclu- 
sions. Yet, they may, and apparently do as we shall see below, indicate a different 
asymptotic behavior, as F — > oo, compared to the simple relations (3.2) and (3.3). 
To clarify this important question we turned to the very limit F = oo. To reach 
this limit we may fix parameters / and /i, and set 7 = u 2 = (see scaling (2.1) and 
(2.2)). In this way we arrive at a conservative system with the simple Hermitian 
Hamiltonian: 

H (ai, CI2) = na\a* 2 + /a* + c.c. = (3.5) 

where the latter equality is due to the initial conditions (2.3). 

A few sample runs of this limit system revealed that all characteristics (2.5) but 
the dimension d of an energy surface now, which is the limit of a chaotic attractor for 
a finite 7, kept growing with time. This suggests an unbounded motion along some 
non-compact energy surface, and moreover with the ever increasing dimensional 
Lyapunov exponent Ai — > 00 for t — > 00. 
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As to the dimension of the energy surface itself it quickly converges to the max- 
imal value d — 4. At the first glance, it appears strange as one would expect this 
surface to be three-dimensional. The explanation of the apparent paradox is the 
following. The metric of the motion M = (+, 0, 0, — ) has now two zero Lyapunov 
exponents: one corresponds to the eigenvector along a trajectoty, as usual, while the 
other one is associated with the eigenvector across the energy surface (the so-called 
marginal instability, see, e.g., Ref.||). It means that this metric characterizes an 
infinitely narrow four-dimentional layer around the energy surface rather than the 
surface itself. Notice that the invariant measure of ergodic motion (the so-called 
microcanonical distribution) on energy surface is determined by the phase-space 
volume in that layer, and not by area on the surface. 

Increasing Ai in the limit (3.5) is at variance with the asymptotics of Eq.(3.2) 
for 7 — > 0. Indeed 

A! = A l7 - M.jUE^j _> , 7-0 (3.6) 

To cope with this difficulty we need a power-law argument in the right-hand side 
of Eq.(3.3) rather than a logarithmic one, F L . To this end, we simply added to the 
latter a power-law term in the form 

F L - F L + exp((3(F L - Fx)) (3.7) 

with two parameters, j3 and Fi. The latter characterizes the crossover to the asymp- 
totic behavior where the Lyapunov exponent is described by the relation 

Fo 7 2/3_1 

The limit behavior in system (3.5) (Ai growth) implies f3 > 1/2. 

There is an interesting possibility to calculate more accurate value of j3 using 
our preliminary empirical data for the limit (3.5). The point is that conservative 
system (3.5) possesses only one quantity of the frequency dimension, the inverse 
time. So, one may conjecture that the variable 1/t plays in conservative system 
(3.5) a role similar to that of 7 in our main dissipative model (1.1). According to 
our preliminary data all the motion characteristics (2.5) but d in the former case 
are growing with time as a power law. Particularly, Xi oc t 1 ^ 3 . By comparison with 
Eq.(3.8) we obtain (3 = 2/3 to be used below. 

With all these improvements we show in Fig. 3 our final fitting the empirical data 
for the attractor dimension and Lyapunov exponents. 

The most interesting result of our studies is finding eventually the two quite 
different chaotic structures with respect to their dependence on the principal pa- 
rameter F. The most of our empirical data belong to the structure with the slow, 
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(3.8) 




10 12 14 16 18 



ln(l+(F-F cr )) 



Figure 3: Same as in Fig. 2 with two more identified attractors (see the Table), and 
the statistical estimate (3.4) for the rest. Two curves are old empirical relations (3.3) 
and (3.2) with a new term (3.7) and fitting parameters: F s = 0.8, F = 7.5, F± = 
11.5. Circles show new data for Ai (see Section 4). 



logarithmic, dependence. It comprises the region from 5F = (SF)i ~ 1 up to the 
crossover 

(5F) 2 « e Fl « 10 5 (3.9) 

Beyond (8F) 2 the chaotic structure becomes 'turbulent' with a relatively fast, power- 
law, dependence on F. The crossover between the two regions is clearly seen in Fig. 3 
from the A 1 (Fj / ) dependence, and can be noticed in cI(Fl) variation as well. So far, 
we have no explanation as to a large value of this crossover. However, notice that 
in the logarithmic scale F L the crossover (F L ) 2 ~ F 1 ~ 10 is not that big. 

Since the logarithmic dependence in Eq.(3.3) becomes linear for small 5F there 
is actually the third region bounded from above by one more crossover (SF) 1 ~ 1. 
This region of threshold chaos is most clearly seen in Fig.l. 
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All the regions can be descibed by unified and relatively simple empirical relations 
(3.2), (3.3), and (3.7) to a reasonable accuracy. An interesting peculiarity of these 
relations is the discontinuity of the first derivative s = dAi(Fi) / 'dFi at d = 3. The 
ratio of the slopes above to below this point 




(3.10) 



is independent of the function itself, similar to Eq.(3.2). 




10 12 14 16 18 

ln(l+(F-F cr )) 



Figure 4: Same as in Fig. 3 for the size R (2.5c) (upper points) and shape S (2.5d) 
(lower points). 

Finally, we display in Fig. 4 the behavior of two geometrical characteristics of 
chaotic attractor, the size R (2.5c) and shape S (2.5d). We did not attempt to 
construct the empirical relations for these characteristics. However, a clear qual- 
itative behavior of those is also interesting and important. First, both R and S 
grow with Fl which is also in agreement with the limit behavior in system (3.5). In 
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other words, the attractor becomes more and more large and elongated. Another 
interesting point is in that the geometrical characteristics are multivalued that is 
specific, at the same value of Fl, for a particular attractor. This confirms that we, 
indeed, have found some different attractors as summarized in the Table. 

4 Conclusion 

In the present work, which continues our previous studies M, we investigate the 
structure of chaotic atractors in the Khalatnikov - Kroyter model of two freedoms 
(1.1). As was shown in Ref. [0, [TJ this simple model describes the birth of the wave 
turbulence (dynamical chaos) in a real very complicated physical phenomenon - 
the dynamics of the second sound in helium II driven by an external force /. The 
model is applicable for weak force only which proved to be sufficient for producing 
chaos. It means that under particular conditions such a few-freedom model can 
represent some features of the many-freedom turbulent structure. The main goal of 
the present studies was to find out what would be peculiarities of the few-freedom 
'turbulence' in a wide range of the parameters, particularly of the driving force. To 
put it other way, we were interested in the structure of the few-freedom chaos as 
compared to that of the multi-freedom one. 

To this end, we have chosen a number of dimensionless characteristics of the 
chaotic structure (2.5), und studied their dependence on a single dimensionless pa- 
rameter F (2.1) in a series of numerical experiments. We had expected some qual- 
itative differences between the few- and multi-freedom chaos, yet the surprise was 
in the nature of the difference which proved to be the dependence of chaos char- 
acteristics on the principal parameter F. While in the multi-freedom chaos this 
dependence is fast (a power law) in a few-freedom chaos it is very slow (logarith- 
mic). Even though our studies were restricted to a particular model we conjecture 
that this principal difference is of a generic nature. Apparently, it is related to the 
strict limitation of the attractor dimension in a few-freedom system (to the interval 
2 < d < 4 in our model). 

In spite of satisfaction with this finding we went on in numerical experiments 
trying to reach as large F as possible in a reasonable computation time (days), and 
it was not in vain! First, we tried to find some indication, if any, of deviations from 
the original empirical relations which could be interpreted as a lower asymptotic 
dimension (d < 4) than assumed above (d = 4). Eventually, we have found the 
deviation but of the opposite sign (see Fig.l and 3)! Even though this is just 
three points we put forward a different explanation - the transition to a new, power 
law, dependence of the attractor characteristics which is similar to that in a multi- 
freedom system (Section 3). 

Currently, we distinguish, altogether, the three regions in the principal parameter 
F. In terms of Ai they are as follows (see Figs.l and 3, and Eqs.(3.2),(3.3),(3.7), 
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and(3.9)): 

(i) the threshold domain 

Ai = -A 3 (d - 2) « i^i . ( F _ Fcr ) ; F cr < F < F cr + 1 (4.1a) 

with a 'turbulent' behavior (cf. Ref. [l|]); 

(ii) the logarithmic region 

F cr + 1 < F < {5F) 2 sa 10 5 (4.16) 

this is our main result revealing a qualitative difference between the few- and multi- 
freedom chaos, and 
(hi) the asymptotics 

Al « • e- 2 ^/ 3 ■ F 2 / 3 , F « oF > (oF) 2 (4.1c) 

-TO 

where the dynamics is similar again to a multi-freedom chaos ('turbulence', cf. the 
first region (4.1a)) provided unbounded motion in the conservative limit (3.5). 

Remarkably, all the three regions can be described by a unified empirical relation 
which is rather simple and reasonably accurate. Particularly, this confirms the 
assumption that (rather than answer to the question if) the dimension of chaotic 
attractor in the model considered does reach the upper bound d — 4 (see Eq.(3.3) 
and above). 

At the time of publication Ref.[0, some preliminary numerical experiments at 
very large F ~ 10 7 (/ ~ 3 x 10 7 ) have presented the first evidence for the transition 
to a power-law dependence in our few-freedom model. It was just three points for 
Ai near the crossover (3.9) marked in Fig.3 by squares. By now, we have collected 
much more empirical data at the crossover and beyond, up to F pa 7 x 10 7 . These 
are also shown in Fig.3 (circles), and perfectly fit, to a less accuracy though, the 
previously found empirical dependence K\{F — F cr ). Each point in the new region 
costed about a full day of work for one processor on the ALPHA supercomputer. 

Unlike this surprising agreement in extrapolation of Ai, the new data for A3 
show a systematic deviation down to < — A 3 >^ 0.75 above crossover (3.9). At the 
moment it is not clear if this decrease is caused by insufficient computation accuracy, 
we had to permit in order to cope with the enormous computation time, or a real 
transition for A 3 is observed, indeed. In any case, the corresponding change in A x 
is negligible in this region (see Eq.(3.2)). 

There is an interesting possibility to compare the unexpected and unexplained 
invariance of A 3 (F) in our model with the data in Ref. || for a close model of two 
freedoms also. In the latter work the fractal geometry of a chaotic attractor was 
studied in detail. For the main set of model's parameters the Lyapunov exponents 
are 

A n = 3.34, 0, -1.79, -5.55 
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with the same metric M (2.4) as in our model, the dissipativity Ar = —4, and the 
attractor dimension d = 3.28 which is well within the logarithmic region (see Fig.3). 
Even more interesting is the ratio 

^ = 2.23 

which differs from our result by 10% only. 

In the very conclusion, we would like to mention that an important question 
how univeral the described structure of the few-freedom chaos might be, remains 
open and certainly deserves further studies. In particular, the above condition for 
a power-law asymptotics (unbounded motion in conservative limit (3.5)) is hardly 
generic. Rather, it is characteristic for a more narrow class of nonlinear models as 
compared to those with the wide logarithmic region. 

Acknowledgments. We are grateful to M. Kroyter and V.V. Parkhomchuk for 
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